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Abstract 

The nonlinear supratransmission is the property of a nonlinear system 
possessing a natural forbidden band gap to transmit energy of a signal 
with a frequency in the gap by means of generation of nonlinear modes 
(gap solitons). This process is shown to result from a generic instability 
of the evanescent wave profile generated in a nonlinear medium by the 
incident signal. 

1 Introduction 

A plane wave scattering onto a medium which possess a natural forbidden band 
gap, is totally reflected if its frequency falls in the gap. Such may not be the case 
if the medium is nonlinear, as the sine-Gordon model considered in ^| where 
the incident wave was shown to generate, above some amplitude threshold, a 
sequence of nonlinear modes propagating in the medium. 

This threshold was given an analytic expression in 2 (still in the sine-Gordon 
case) and the property of the nonlinear medium to transmit energy in the gap 
has been called nonlinear supratransmission. It has been shown then to hold 
for different models (nonlinear Klein-Gordon, double sine-Gordon, Josephson 
transmission lines), and has been experimentally realized on a chain of coupled 
pendula |3J. 

Energy transmission by gap soliton generation is a phenomenon that has 
been experimentally discovered in Bragg media with Kerr nonlinearity 4 . In 
that case the underlying simplified model is the coupled mode system jS] and 
it has been recently demonstrated that the very mechanism of the switch from 
a Bragg mirror to a medium with high transmissivity is the nonlinear supra- 
transmission [jy. Again there the theory allows for an analytic expression of 
the threshold amplitude, based on the expression of the soliton solution of the 
coupled mode system found in [7]. 

We demonstrate here that the process of nonlinear supratransmission is a 
natural consequence of a generic instability of the evanescent profile generated 
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in a nonlinear medium when it is constrained at one end by a periodic boundary 
driving of frequency in a forbidden band gap. 

2 General context 

Let us consider the Lagrangian of a scalar field u(x, t) in the half-plane (x, t) € 

\u 2 t - ^ul - <f>(u) . (1) 



4>'(u) = , (2) 

where prime denotes differentiation with respect to the argument. The wave 
equation @ is submitted to Dirichlet boundary conditions 

u(0, *)=/(<), «((»,*) = 0, (3) 

associated with Cauchy initial data u(x, 0) and ut(x, 0), which can be vanishing 
or adapted to the boundary forcing. Note that this approach can be extended 
to Neumann boundary conditions. 

From now on we adopt the quite generic assumtion of a potential energy <p(u) 
which possess a local minimum locally symmetric. By a convenient choice of 
the reference frame, this minimum can be taken to be located at (</>, u) = (0, 0). 
Then the symmetric requirement means that the potential is locally an even 
function of u. Moreover, by a convenient choice of the units of x and t we can 
always scale the value of (f>"(0) to the value 1. In summary the potential <p(u) 
obeys the following requirements 

0(0) =0, ^'(0)=0, 0"(O) = 1, <T(0) = 0. (4) 

Note that the local parity assumption (f>"'(0) — is not essential but simplifies 
the expressions. 

As a consequence the equation @, by Taylor expansion of </>'(u) about the 
minimum u = 0, has the linear limit 

utt - u xx + u = , (5) 

with dispersion relation 

uj 2 = 1 + k 2 , (6) 

having the gap [—1, +1] in the frequency. 

The second restiction is on the class of boundary datum fit) which we take 
to be a periodic forcing with a frequency belonging to the forbidden band gap. 
More precisely we assume 

u(0,t) = A cos(fii) , |0| < 1 . (7) 



/■oc 

C = / dx Liu,ut,u x ) , L = 
Jo 

The equation of motion 5C = reads 

x > : u tt - Uxx + 
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This problem is quite generic and describes e.g. the scattering of a monochro- 
matic wave into a medium in one of its forbidden band gaps. 

For the linear limit JSJ, the boundary condition (7J) would generate the so- 
lution (x > 0) 

u (x, t)=A cos(ftt) e~ Kx , fl 2 = 1 - K 2 , (8) 

called the linear evanescent wave. This solution is a reasonable approximate 
solution of the nonlinear evolution as long as the amplitude A is small enough. 

We show hereafter that in the nonlinear case, under some additional simple 
conditions on the potential 4>{u), the medium can present nonlinear supratrans- 
mission as a result of a fundamental instability generated by the evanescent 
wave, above some threshold for the amplitude A. 



3 Perturbation scheme 



Let us start with the simplest nonlinear correction resulting from the assumption 
Q, namely the nonlinear wave equation 



1 o 

u tt ~ u xx + u + -au = , 
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where the nonlinear factor a is defined by 

_ d 4 cj> 



(9) 



(10) 



u=0 
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Figure 1: Plot of the actual nonlinear evanescent profile ui(x) (full line) com- 
pared to the linear one uq(x) (dashed line) in the static case f2 = for A = 3.4. 
The dots are the results of a numerical simulation of 

It is essential to remark at this step that the first effect of the nonlinearity 
on a boundary driving (JJJ) is to deform the (linear) evanescent profile JHJ) . Such 
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a deformation can be illustrated analytically and numerically in the static case 
Q = in J7J). Indeed the equation with a = — 1 for u(x) submitted to 
u(0) = A, possess the solution 

= u7~~ x . ^0 = -|arccosh(^^)| . (11) 

cosh(ir — Xq) A 

This is the profile displayed on the figure ^ with A = 3.4 (full line) in perfect 
agreement to numerical simulations (dots). The numerical method is explained 
in the appendix and here we have solved with a boundary datum u(0, t) 
slowly increasing from to the value A = 3.4 and the solution (represented 
by the crosses in fig. is plotted long after the boundary datum has been 
setteled. We have plotted also for comparaison to the exponential (linear) shape 
(|2J namely uq(x) — Ae~ x (dashed line). 

In the dynamical case (fi ^ 0) we denote by ui(x, t) the deformed evanescent 
profile, solution of 10 , hence required to obey the boundary value Q , namely 

ui(0,t) = A cos(ffi) . (12) 

Although this solution u\{x, t) may not be known analytically, we shall actually 
discover that the worst approximation of it, namely the linear evanescent profile 
uq(x, t) itself, still generates the seeked instability, hence providing a qualitative 
understanding of the generation of nonlinear supratransmission. 



4 Instability criterion 

The second effect of the boundary driving is obtained by studying the pertur- 
bation of the exact solution ui(x, t) of © as 

u(x, t) — ui (x, t) + erj(x,t) , e <C 1 . (13) 

At order e, the evolution for the perturbation n then reads 

z>0, rht-Vxx + rill + V] =0 , V(x,t) = -au\{x,t) . (14) 

This linear wave equation in the potential V on the semi-axis x > is our main 
tool and contains all the ingredients required for a qualitative description of the 
onset of instability. 

The linear equation l|14|) can be thought of as a two-dimensional scattering 
problem on the half-plane for a potential which is periodic in the infinite t- 
dimension and localized around x = in the positive cc-direction. The study of 
this mixed periodic-localized 2D scattering problem is a tough technical question 
that we report in a future investigation. 

However we show hereafter that the instability can be qualitatively under- 
stood by considering the static case where V is assumed to depend only on the 
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space variable x. The scattering problem can then be reduced to the radial 
S-wave Schrodinger spectral problem 

fe(C, x) + [(A 2 - 1) - V(x)] V(C, x)=0, (15) 

for the spectral parameter £ and cigcnfunction ip(C, x) defined by 

C 2 = A 2 -1, r](x,t) =e iXt ^x) . (16) 

We argue then that equation (|14fl develops an instability if the potential V 
possess a bound state, say £ = ip, with sufficient energy, namely if 

C 2 = -p\ v > 1 => A 2 < , (17) 

which causes the first order correction 77(2;, f) to diverge in time. In particular 
it is clear from (|14|l that the instability will develop only when a < as a 
necessary (not sufficient) condition to have a bound state for V. In order to 
allow for a comparaison with sine-Gordon, and because the value of a in JSJ 
can be scaled with A, we set 

a=-l. (18) 

The potential V(x) is obtained from the deformed evanescent profile U\(x, f) 
by considering its extremum on a period, namely we take 

V(i) = -^maxK(i,t)}. (19) 

We now illustrate the obtained expression for the threshold of instability for 
two particular choices for u\, first the linear evanescent profile (JHJ, second a 
standard sech-type profile. Both cases allow for an explicit computation of this 
threshold and we compare in figure these results to numerical simulation of 
the PDE ©. 



5 The exponential potential 



The choice Aexp[—Kx) for max{u 2 (ir, t)} provides from l|15(l the following spec- 
tral problem on the half-line 



x > , ip xx ((,x) + 



1 



C 2 + -A z e 



2-IK'jl 



lK£a;)=0 



(20) 



The two (positve, real-valued) parameters are the amplitude A and the wave 
number K = yl — TP. 

The direct scattering problem [HI E] essentially consists in finding the Jost 
solution of (|20|l which can be written 



f(C,x)=e«* Fi(l-iL 



A 2 



,-2K 



K' 8K 2 



(21) 



5 




Figure 2: Comparaison of the threshold of instability for the full equation 
(crosses) with the theory obtained from 11411 with a = — 1 for u\ = Acxp[— Kx] 
(full line left figure) and u\ = A/ cosh(i^x) (full line right figure). 



This solution is usually written in terms of the Bessel function of the first kind of 
order —iC^jK. We find it more convenient to express it in terms of the generalized 
hypergeometric function q-Fi ■ Indeed it readily provides the asymptotic behavior 

f(C,x)e-** — 1 , (22) 

x— >oo 

and the Jost function 

/(CO) =0^(1-^;-^) . (23) 

The bound states are the zeroes of this Jost function in the upper complex em- 
plane, actually located on the imaginary axis as the potential is real valued. Our 
instability criterion l|17|l then tells that an instability for rj(x, t) will develops as 
soon as the amplitude A reaches the value A s for which the first bound state 
reaches the value C = i- This occurs for each value of the parameter K and thus 
the threshold A S (K) is defined by the implicit equation 

oW + i;-^)=0. (24) 

We have solved this equation numerically and the result is plotted in terms of 
= y/i - K 2 in the left figure 

It is compared there with the results of the numerical solution (see ap- 
pendix) of © (with a = — 1) submitted to the boundary driving (JJJ, namely 
u(0,t) = A cos fit. Although the theoretical threshold does not fit perfectly the 
simulations, the principle of the nonlinear supratransmission is qualitatively un- 
derstood. The discrepancy results both from the asumption of an exponential 
profile for ui(x,t) and a static approximation of it. 
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6 The Eckart potential 

A better fit is obtained by assuming instead of the exponential profile the follow- 
ing sech-type profile obtained from the solution (|if() at its maximum amplitude 
(i.e. for xq — where the instability is expected) 

Ul = ^> K = VT^. (25) 

The presence of a periodic boundary datum is taken into account by the param- 
eter K (wave number). This choice is expected to give good results for small ft 
and it leads to the spectral problem 



x > , ip xx ((,x) + 



c 2 + ^- 



f A 2 



i/j(C,x) = 0. (26) 



2 cosh 2 Kx_ 

This potential belongs to the Eckart class and its Jost solution reads 

f(C, x) = e< x (1 + e- 2Kx f 2 F 1 (/3, /3 ~ ^ 1 " ijf, ~^ Kx ) > (27) 
with the parameter 

P = I f 1 - 71 + 2A 2 /K 2 ) . (28) 



/(C 0) = T^T-fzT^ > (3°) 



2 

This solution does obey the asymptotic (|22|l and the Jost function is 

/(C, 0) = 2^ 2 F 1 (/3, 0-i^;l-ij^; -1) . (29) 

According to the instability criterion (|f 7|l we then seek the value A s for the 
threshold amplitude generating a bound state (zero of this Jost function on the 
upper imaginary £-axis) at the value £ = i. From ^Uj we can write also l|29[) as 

2*/*y?r(i-4) 

\ 2 2K I \ 2 '2if / 

where T(v) is the standard Euler Gamma function that possess poles at v = 
0,-1,.-.. 

The bound states are thus given by the poles of the function r(^i2 — i4jA 
on the positive imaginary £-axis, and the first bound state reaches the point 
£ = i as soon as /3 verifies 

i±^ + ^=0. (31) 
2 2K y ' 

Replacing there the value l|28l) of the parameter j3, the threshold of instability 
A a can then be expressed in terms of f2 = vl — if 2 as 



^ 2 = 4(1 -ft 2 ) +6Vl- ^ 2 + 2 . (32) 

This is the curve represented in the right graph of figure [21 We see then that 
the sech-like profile provides as expected a much better fit of the instability 
threshold, especially for small ft where the static approximation becomes more 
accurate. 
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7 The sine- Gordon case 



So far we have demonstrated the existence of an onset of instability in the 
nonlinear simple model 10 (with a < 0). In that case the instability causes 
the solution to diverge because the related potential 4>{u) = u 2 /2 — u 4 /24 is not 
confining for large u. 

To illustrate the nonlinear supratransmission as an effect of the instability, 
we consider now the sine-Gordon case <f)(u) = 1 — cos(u) treated in 0E], namely 



x > 



Utt — u xx + sin u = , u(0, i) = A cos fit 



(33) 



In this situation, the solution induced by the boundary datum of amplitude 
below the threshold of instability is the non-travelling breather 



Ub{x,t) = 4arctan 



K 



cos(fit) 



fl cosh(K(x — Xb) 



Xb < , 



(34) 



where Xb is the breather center. 

The onset of nonlinear supratransmission is reached, at given frequency fi, 
for the maximum amplitude of Ub(0,t), obtained of course for xt, = 0. Let 
v,2(x,t) denote the solution Ub for Xb = 0, namely 



U2(x,t) = 4arctan 



K cos(fit) 



n cosh(Kx 



(35) 



The first order perturbation ei](x, t) about the solution ui obeys from (|33|l the 
linear evolution 



x>0, f]tt-r) xx + r)[l + V] = , V(x, t) = cos(it 2 ) - 1 



(36) 



that replaces equation (|14|l used precedingly. The potential V can also be 
written 



V(x,t) 



, K 2 cos 2 fit 
' n 2 cosh 2 Kx 



1 



K 2 cos 2 fit 



fl 2 cosh 2 Kx 



(37) 



Our purpose here is to show that the threshold of instability of the linear 
evolution (|36(l with the stationary approximation ()19(l . namely 



V(x) 



D 2 



cosh Kx 



1 



B 2 



1 -2 



cosh Kx 



(38) 



matches reasonably the theoretical prediction, namely the value B s of the pa- 
rameter B given from (|37f) by 



K 



Vi -fi 2 
fi 



(39) 



The drawback here is that the potential l|38ll does not lead to an explicitely 
solvable spectral problem (|15|) . This is why we have used numerical simulations 
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Figure 3: Onset of instability for the parameter B of V(x) given by numerical 
simulations of i|3t)|) with 1|38|) (dots) compared to the prediction (|39|) . 

and shown by figure|3]that the stationary limit does produce a consistent result: 
the threshold of instability of the PDE l|36[l obtained by varying B, at each fixed 
r2, fits suprisingly well the exact prediction (|39"|) . 

Note finally that in the sine-Gordon case we have an analytic expression of 
the threshold amplitude that generates nonlinear supratransmission 2 , namely 

A S (Q) = 4arctan[ v / l - W/tt] . (40) 

This is simply the maximum amplitude of the breather (|35[1 . which is indeed 
reached when B in 1|38|) takes the value 1)39(1 . Above this value the instability 
comes into play and generates a nonlinear mode that propagates in the medium 
carrying energy, what have been called nonlinear supratransmission and have 
been experimentally realized on a chain of pendula 3 . 

8 Conclusion 

We have demonstrated that a nonlinear medium, under some generic condi- 
tions on the structure of the nonlincarity (existence of a forbidden band gap, 
local minimum for the potential, ...) may create an instability of the evanes- 
cent profile generated by a periodic boundary driving at a frequency in the 
gap. This instability is the generating process at the origin of nonlinear supra- 
transmission, when the nonlinear medium supports nonlinear modes (breathers, 
solitons, kinks, ...). 

The instability has been described by using a stationary limit for the evanes- 
cent profile which has allowed us to prove that the occurence of the instability 
is related to the appearence of a bound state with sufficient energy, according 

to jnji. 

We have dealt here with a scalar wave- like equation and it is worth mention- 
ning that the process works also for different types of systems like for instance 
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the coupled mode equations in Bragg media The generating instability in 
this case will be considered in a future work. 

A Numerical simulations 

In order to integrate numerically the partial differential equations to deal with, 
we have used two different approaches depending wether the PDE is linear or 
not. 

In the nonlinear case, such as equation (|2J), we have used an explicit second 
order scheme in space and solved for time evolution as a set of coupled differ- 
ential equations. More precisely we have replaced J2J with the set of ODE's 
(overdot means time differentiation) 

u n - j^(u n +i - 2u„ + Un-i) + 4>'{un) = , n = 0..N , (41) 

where h is the chosen grid step and u n (t) = u(nh,t). 

The boundary datum is smoothly setteled by assuming 

u v(t) = A- [1 + tanh(p(t - to))] cosilt , (42) 

where p and to are the parameters to select. A semi-infinite line is simulated 
by an absorbing end consisting of a damping term 7u„ included in the equation 
where the intensity 7 slowly varies from to 1 over the last 40 cells. 

The results presented in figures and |21 are obtained with N — 100 spatial 
mesh points, h = 0.2 grid spacing over a time of integration t m = 200. The 
parameters for the boundary field are p = 0.2 and to = 20. The resulting set 
of coupled ODEs is then solved through the subroutine dsolve of the MAPLE8 
software package that uses a Fehlberg fourth/fifth order Runge-Kutta method. 

In the linear case, namely for equation (|36l) . we have used the pdsolve 
routine of MAPLE8 over a length L — 10 and a space-time step of 0.1. It uses a 
centered implicit scheme, of second order in space and time, with a first order 
boundary condition. To test the stability of an equation like H36|) we have solved 
it for the following set of initial-boundary values 

?7(0, t) = 0.01 , r](L, t) = 0.01 , r)(x, 0) = 0.01 , T) t (x, 0) = . (43) 

In the stable region we observe small amplitude oscillations that grow when the 
potential parameter get close to the instability threshold. Above this threshold 
the instability manisfests as an exponential growth of the solution. The points 
that result in the figure [21 denote the first value of the parameter B for which 
instability settles in, with an absolute precision of 10 _1 . 

References 

[1] J-G. Caputo, J. Leon, A. Spire, Phys Lett A 283 (2001) 129 



10 



[2] F. Gcnict, J. Leon, Phys Rev Lett 89 (2002) 134102 

[3] F. Geniet, J. Leon, J Phys Cond Mat 15 (2003) 2933 

[4] D. Taverner, N.G.R. Broderick, D.J. Richardson, R.I. Laming, M. Ibsen, 
Opt Lett 23 (1998) 328 

[5] C. Martijn de Sterke, J.E. Sipe, Phys Rev A 38 (1988) 5149 

[6] J. Leon, A. Spire, Gap soliton formation by nonlinear supratransmission 
in Bragg media, preprint Montpellicr 2003 

[7] A.B. Aceves, S. Wabnitz, Phys Lett A 14 (1989) 37 

[8] R.G. Newton, Scattering theory of waves and particles, McGraw-Hill 
(New- York 1966) 

[9] K. Chadan, P.C. Sabatier, Inverse problems in quantum scattering theory, 
Springer- Verlag (Berlin 1989) 

[10] A. Erdelyi, Higher transcendental functions, vol 1, McGraw-Hill (New- 
York 1953) 



11 



